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Abstract 



Using the D0 detector, we have studied events produced in pp collisions 
that contain large forward regions with very little energy deposition ("ra- 
pidity gaps") and concurrent jet production at center-of-mass energies of 
y/s = 630 and 1800 GeV. The fractions of forward and central jet events as- 
sociated with such rapidity gaps are measured and compared to predictions 
from Monte Carlo models. For hard diffractive candidate events, we use the 
calorimeter to extract the fractional momentum loss of the scattered protons. 
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Typeset using REVTeX 



Inelastic diffractive collisions are responsible for 10-15% of the pp total cross section and 
have been described by Regge theory through the exchange of a pomeron Diffractive 
events are characterized by the absence of significant hadronic particle activity over a large 
region of rapidity or pseudorapidity (rj = — ln[tan(|)], where 9 is the polar angle relative to 
the beam). This empty region is called a rapidity gap and can be used as an experimental 
signature for diffraction. Recent interest in diffraction has centered on the possible partonic 
nature of the pomeron in the framework of quantum chromodynamics (QCD), as suggested 
by Ingelman and Schlein ||. Hard single diffraction (HSD), which combines diffraction and 
a hard scatter (such as jet or H^-boson production), can be used to study the properties of 
the pomeron. 

The partonic nature of the pomeron was first inferred by the UA8 experiment |§ at the 
CERN SppS collider at ^/s = 630 GeV from studies of diffractive jet events. Recent analyses 
of diffractive jet production |§-|| and diffractive W^-boson production @ are consistent with 
a predominantly hard gluonic pomeron, but measured rates at the Fermilab Tevatron are 
several times lower than predictions based on data from the DESY ep collider HERA ||. 
In this Letter we present new measurements of the characteristics of diffractive jet events, 
and of the fraction of central and forward jet events that contain forward rapidity gaps 
("gap fraction") at center-of-mass energies ^/s = 630 and 1800 GeV. These measurements 
augment previous results from the CDF collaboration on the gap fraction for forward jets 
at y/s = 1800 GeV (3| and place further constraints on diffractive models. 

In the D0 detector ||, jets are measured using the uranium/liquid-argon calorimeters 
with an electromagnetic section extending to \r]\ <4.1 and coverage for hadrons to \r]\ <5.2. 
Jets are reconstructed using a fixed-cone algorithm with radius 1Z = y/Arj 2 + A0 2 = 0.7 (0 
is the azimuthal angle). The jets are corrected using standard D0 routines for jet-energy 
scale |TIJ, except that there is no subtraction of energy from spectator parton interactions, 
since these are unlikely for diffractive events. 

To identify rapidity gaps, we measure the number of tiles containing a signal in the L0 
forward scintillator arrays {n^), and towers (Arj x A<p = 0.1 x 0.1) above threshold in the 
calorimeters (t^cal)- The L0 arrays provide partial coverage in the region 2.3 < \r]\ < 4.3. A 
portion of the two forward calorimeters (3.0 < \r]\ < 5.2) is used to measure the calorimeter 
multiplicity, with a particle tagged by the deposition of more than 150 (500) MeV of energy 
in an electromagnetic (hadronic) calorimeter tower. The thresholds are set to give negligible 
noise from uranium decays, while maximizing sensitivity to energetic particles |TT . 



For i/i = 630 and 1800 GeV, we use triggers which required at least two jets with 
transverse energy Et > 12 or 15 GeV (see Table |) to study the dependence of the gap 
fraction on jet location. The forward jet triggers required the two leading jets to both have 
i] > 1.6 (or rj < —1.6), while the central jet triggers had an offline requirement of \r]\ < 1.0. 
These data were obtained during special low luminosity runs, with typical instantaneous 
luminosities much less than 1 x 10 30 cm' 2 s _1 . At each y/s, we also implemented the so- 
called single veto trigger (SV), a dijet trigger that required a rapidity gap on one side (using 
the L0 detector). The SV trigger was used to obtain large samples of single diffractive 
candidate events. The events in the final data samples all have a single pp interaction 
requirement, a vertex position within 50 cm of the center of the interaction region, and two 
leading jets that satisfy standard quality criteria [Q. The number of events in each of the 
final data samples and the integrated luminosities (£) are given in Table Q. 
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TABLE I. Attributes of the final data samples. 



Data Sample Jet |ry| Jet -£/r(GeV) C (nb 1 ) Events 

1800 GeV Forward > 1.6 >1~2 62^9 50852 

1800 GeV Central < 1.0 > 15 4.55 16567 

630 GeV Forward > 1.6 > 12 16.9 28421 

630 GeV Central < 1.0 > 12 8.06 48123 

1800 GeV SV - > 15 5700 170393 

630 GeV SV - > 12 529 64772 



The nL0 versus ncAL distributions for central and forward jet events at y/s = 630 and 
1800 GeV are shown in Fig. [I]. For forward jet events, these quantities are defined by the rj 
region on the side opposite the two leading jets, while for central jet events they are defined 
by the forward r\ interval that has the lower multiplicity. The distributions display a peak 
at zero multiplicity (n C AL = ^L0 = 0), in qualitative agreement with expectations for a 
diffractive component in the data. 
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FIG. 1. Multiplicity distributions at y/s = 1800 GeV for (a) forward and (b) central jet events, 
and at y/s = 630 GeV for (c) forward and (d) central jet events. 

The gap fraction is extracted from a two-dimensional fit to the lego plot of versus 
^-cal- The non-diffractive (high multiplicity) background is fitted in the signal region using 
a four-parameter polynomial, and the signal is fitted with a falling exponential, as suggested 
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by Monte Carlo [|Tl|]. Figure ^| shows the multiplicity distribution from Fig. |](a), an d 
the resulting fitted signal, fitted background, and normalized distribution of pulls ( [data- 
fit] /a/AT). All distributions have adequate fits, with x 2 /dof < 1.2. 

Table [D] shows the gap fractions obtained for the four event samples. The values range 
from (0.22 ± 0.05)% for central jets at y/7s = 1800 GeV, to (1.19 ± 0.08)% for forward jets at 
y/s = 630 GeV. Uncertainties are dominated by those on the fit parameters. Additional small 
uncertainties from the dependence on the range of multiplicity used in the fits were added 
in quadrature. Potential sources of systematic error, such as the number of fit parameters, 
jet energy scale, trigger turn-on, tower threshold, luminosity, residual noise, and jet quality, 
yield only negligible variations in the gap fractions |p"T| . 

Table || shows that the gap fractions at ^/s = 630 GeV are larger than gap fractions 
at yfs = 1800 GeV and that gap fractions for forward jets are larger than for central jets. 
Table |T| also lists predicted gap fractions for several possible pomeron structure functions 
(discussed below). 
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FIG. 2. The (a) data from Fig. ||(a), and corresponding (b) fitted signal, (c) fitted background, 
and (d) normalized pull distributions. 

We compare the data to Monte Carlo (MC) simulations using the hard diffractive event 
generator POMPYT |]13" |, which is based on the non-diffractive PYTHIA |L4|] program. In 
POMPYT, a pomeron is emitted from the proton with a certain probability (called the flux 
factor [0]), and has a structure functions s(/3), where j3 is the fractional momentum of 
the pomeron carried by the hard parton. We used the standard Donnachie-Lanshoff flux 
factor [[DJ in this analysis and compare our data to four structure functions: (i) "hard gluon," 
a pomeron consisting of two gluons, s(/3) oc (3{1 — /?); (ii) "flat gluon," s(/3) oc constant; 
(iii) "soft gluon," s{(3) oc (1 — /?) 5 ; and (iv) "quark," the two-quark analog of (i). In each 
case, the gap fraction is defined as the cross section for jet events with a rapidity gap based 
on POMPYT divided by the jet cross section from PYTHIA Many uncertainties, such as the 
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choice of proton parton densities, cancel in the ratio. An MC version of the fitting method 
is applied to correct for diffractive events that fail the gap selection criteria. By applying 
the appropriate correction factor (which ranges from a few per cent for soft gluon central 
jets to about 80% for hard gluon forward jets) to each MC prediction and comparing to the 



data, we make no assumptions about which model (if any) is correct ||T6 |. 

Monte Carlo gap fractions are shown in Table [0]. The systematic uncertainties are 
typically dominated by the difference in energy scale between data and Monte Carlo, but 
also include uncertainties from the fitting procedure. We observe that rates for harder gluon 
structures are far higher than supported by data, while the quark structure is in reasonable 
agreement with the data. The quark structure, however, has previously been shown to 
predict an excessive rate of diffractive IF-Bosons 0. 

A hard gluonic pomeron is capable of describing previous measurements pH7j, if com- 



bined with a flux factor that decreases with increasing yfs [|TTf| . The ratios of gap fractions 
shown in the lower half of Table [H] provide new information, since the flux factor cancels for 
the same y/s, and dependence on the flux factor is reduced for different y/s. The ratios for 
jets with \rj\ > 1.6 to jets with < 1.0 show clear disagreement between the data and pre- 
dictions for a hard-gluon pomeron structure, despite this cancellation. A gluon-dominated 
pomeron containing both soft and hard components, combined with a reduced flux factor, 
could describe all the data samples. 

TABLE II. The measured and predicted gap fractions and their ratios. 



Sample 



Data 



Gap Fractions 
Hard Gluon Flat Gluon 



Soft Gluon 



Quark 



1800 GeV \rj\ > 1.6 
1800 GeV \rj\ < 1.0 
630 GeV \rj\ > 1.6 
630 GeV Inl < 1.0 



(0.65 ±0.04)% (2.2 ±0.3)% (2.2 ± 0.3)% 

(0.22 ±0.05% (2.5 ±0.4)% (3.5 ± 0.5)% 

(1.19 ±0.08)% (3.9 ±0.9)% (3.1 ±0.8)% 

(0.90 ±0.06)% (5.2 ±0.7)% (6.3 ±0.9)% 



(1.4 ±0.^70 

(0.05 ±0.01)% 

(1.9 ±0.4)% 
(0.14 ±0.04)% 



(0.79 ±0.12)% 
(0.49 ± 0.06)% 
(2.2 ± 0.5)% 
(1.6 ±0.2)% 



Ratios of Gap Fractions 



630/1800 H > 1.6 1.8 ±0.2 1.7 ±0.4 1.4 ±0.3 1.4 ±0.3 2.7 ±0.6 

630/1800 M < 1.0 4.1 ±0.9 2.1 ±0.4 1.8 ±0.3 3.1 ± 1.1 3.2 ± 0.5 

1800 | r? I > 1-6/M < 1-0 3.0 ±0.7 0.88 ±0.18 0.64 ±0.12 30. ± 8. 1.6 ±0.3 

630 \rj\ > 1.6/|r/| < 1.0 1.3 ±0.1 0.75 ±0.16 0.48 ±0.12 13. ± 4. 1.4 ±0.3 



The characteristics of the HSD events were examined using the high statistics SV trigger. 
We plot in Fig. || the distributions of the number of jets, the E T - weighted rms jet widths, the 
A<f) between the two leading jets, and the relative ratio of diffractive to non-diffractive events 
as a function of the average Et of the two leading jets, for central jets at \fs = 1800 GeV. 
The solid lines in Fig. [^(a)-(c) correspond to the distributions for HSD candidate events 
(^cal = n L& = 0), and the dashed lines show the distributions for non-diffractive events 
(^cal > and > 0). These plots show that the diffractive events appear to have less 
overall radiation. Figure |3](d) indicates that there is little dependence of the gap fraction on 
average jet Et- The MC samples (not shown) have characteristics similar to the data. 
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FIG. 3. Distributions of the (a) number of jets, (b) jet width, (c) A0 between leading jets, 
for central diffractive (solid) and non-diffractive (dashed) jet events at y/s = 1800 GeV. (d) The 
relative ratio of diffractive to non-diffractive events as a function of the average Er of the two 
leading jets. 



Finally, we measure the fractional momentum loss of the proton £, denned as [|T^]: 

((U) 

V s i 

where the summation is over all observed particles. The outgoing scattered proton or an- 
tiproton (and the rapidity gap) is defined to be at positive rj. Equation (1) weights heavily 
the well-measured central region near the rapidity gap, while particles that escape down the 
beam pipe at negative rj give a negligible contribution. Using POMPYT events, where £ can 
be determined from the momentum of the scattered proton, we have verified that Eq. (1) is 
reliable at both values of y/s and for different pomeron structures. A scale factor (2.2 ± 0.3) 
derived from Monte Carlo is used to convert £ measured from all particles to that from 
just electromagnetic calorimetric energy depositions [fIT}| . The £ distributions for forward 
and central jets at y/s = 630 and 1800 GeV are displayed in Fig. £|, with the shaded region 
showing the variance in the distribution due to energy scale uncertainties. Energy-scale un- 
certainties result in a shift in £ such that if the true distribution were below the histogram 
at small £, it would be above the histogram at large £. 

The £ distributions show the expected kinematic behavior of diffraction (M = y/^s, 
where M is the mass of the diffractive system), peaking at larger £ for central jets than for 
forward jets. Forward and central jets at y/s = 630 GeV also peak at larger £ values with 
respect to the corresponding distributions at y/s = 1800 GeV, since for fixed diffractive mass, 
smaller y/s implies larger £. Even though pomeron exchange is thought to dominate only 



9 




0.025 0.05 0.075 0.1 0.025 0.05 0.075 0.1 

(a) % (b) k 




0.05 0.1 0.15 0.2 0.05 0.1 0.15 0.2 

(c) k (d) k 

FIG. 4. The £ distributions for yfs = 1800 GeV (a) forward and (b) central jets and for 
yfs = 630 GeV (c) forward and (d) central jets, using the SV trigger with uqal = "-L0 = 0. 
The shaded region shows the variance in the distribution due to energy scale uncertainties (see 
text). 

for £ < 0.05, the trends of the £ distributions can be reproduced by POMPYT. Without the 
observation of the scattered proton, the interpretation of these large £ rapidity gap events 
is uncertain. 

We have measured properties of hard single diffraction at \fs = 630 and 1800 GeV with 
jets at forward and central rapidities. The gap fractions have been measured without ap- 
plying model-dependent corrections. Within the Ingelman-Schlein model, our data can be 
reasonably described by a pomeron composed dominantly of quarks. For the model to de- 
scribe our data as well as previous measurements, a reduced flux factor convoluted with a 
gluonic pomeron containing significant soft and hard components is required. We have also 
measured the fractional momentum lost by the scattered proton and found it greater than 
typically expected for pomeron exchange. 
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